Mean— field results on the Anderson impurity model out of equilibrium 
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We investigate the mean-field phase diagram of the Anderson impurity model out of equilibrium. 
Generalising the unrestricted Hartree-Fock approach to the non-equilibrium situation we derive and 
analyse the system of equations defining the critical surface separating the magnetic regime from the 
non-magnetic one. An exact analytic solution for the phase boundary as a function of the applied 
voltage is found in the symmetric case. Surprisingly, we find that as soon as there is an asymmetry, 
even small, between the contacts, no finite voltage is able to destroy the magnetic regime which 
persists at arbitrary high voltages. 
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The field of 'Kondo physics' sprung from the crucial 
insight by Anderson into the formation of local magnetic 
moments in metals Q. In Ref.p], Anderson formulated 
a very simple model, now known as Anderson's impu- 
rity model (AIM), and showed that there is a critical 
value of Coulomb coupling (U c ) above which a local mag- 
netic moment spontaneously forms (in the framework of 
a mean- field approximation). It was understood later 
(see 0, especially Chapter 5, and references therein), by 
using more sophisticated methods, that U c is in fact a 
crossover rather than a sharp 'phase transition'. Never- 
theless, the Anderson's original work is a milestone 
for our understanding of the 'standard' Kondo problem. 
It still serves as a foundation for various more refined 
techniques, a notable example being the 'local moment 
approach' (LMA)|jJ. 

In recent years, considerable theoretical interest in the 
non-equilibrium AIM developed, fuelled by experimen- 
tal results on non-equilibrium transport in quantum dots 
where Kondo effect signatures were observed ii. The 
low bias voltage (V) limit of the AIM has been almost 
completely understood [fj, |j] (perhaps the precise mech- 
anism of 'dephasing' still being in need of some clarifica- 
tion ||)- When the voltage becomes larger than some 
characteristic scale (like the Kondo temperature), the 
situation is less clear. While there are some interest- 
ing results in the literature jjj, no consistent theoretical 
method seems to be available to tackle the problem be- 
yond the perturbation expansion in the on-site interac- 
tion, see e.g. Ref.|lll|. In this paper we provide a start- 
ing point for one. Our objective is to generalise Ander- 
son's original analysis [unrestricted Hartree-Fock (HF) 
approach] to the out of equilibrium model. Note that 
our starting point is complementary to that suggested 
in Ref.|lO|. which is appropriate for an expansion in the 
electron tunnelling amplitudes. 

We start with the standard AIM Hamiltonian. It con- 
tains a single site with energy A described by electron 
creation and annihilation operators (f a , d a , where a de- 
notes the spin orientation. This site is coupled to two 
non-interacting electron systems i = 1,2, which we shall 



also call 'leads'. Their dynamics is governed by the 
Hamiltonians Hosier]) where ipi a are the corresponding 
electron field operators. The coupling of the leads to 
the impurity site takes place via hopping with energy- 
independent amplitudes 7r,2, 



cr=T,i i=l,2 

+ 744^(0) + Vl(0)da]}} + Ud\d^ d\d t . (1) 



The last term describes the on-site repulsion on the impu- 
rity. The cornerstone of the HF approach is the following 
approximation for this last term, 



U(n^d[d l +n i djd T ) 



(2) 



where n = (d\d a ) are the level population probabili- 
ties which have to be found self-consistently. The self- 
consistency equations are basically the population prob- 
abilities as functions of each other. In order to calculate 
them we first fix both n a . Then the full Hamiltonian sep- 
arates with respect to spin degrees of freedom and can 
easily be diagonalised keeping in mind that n a are pa- 
rameters. Since we are dealing with a non-equilibrium 
situation, the calculation of the population probabilities 
can only be accomplished using suitable modifications of 
the conventional equilibrium techniques. We decided to 
use the Keldysh formalism |l2|,[l3|. Let us first define the 
local (taken at the site adjacent to the impurity) Green's 
function of the lead electrons (from now on we set the 
spacial coordinate to x = and drop it hereafter), 



G ia {t-t') = -i{T c ^{t)4At')) 



(3) 



where Tc denotes the time ordering operation along the 
closed Keldysh contour C consisting of the forward path 
CL, running from — oo to oo, and the backward path 
going the other way round. After an appropriate 
placement of the time variables, one can recover all usual 
Keldysh Green's functions from the generalised function 
(J3J, we use notation of Ref.0. We shall use only 
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the bare Green's functions of the lead electrons, i. e. 
those calculated without the tunnelling ji = [the up- 
per (lower) sign and subscript 1(2) corresponds to the left 
(right) lead], 



. Po 



i( 2) » = G^ ) » = -i^Bgn( W ±V/2) 



1(2)<7 

+- 

1(2)<7 



-ip e{uj±V/2) 



(4) 



where po is the constant density of states in the leads and 
6 denotes the Heaviside step function. 

Green's functions of the impurity level are defined anal- 
ogously and are denoted by D. The function D~ + (t — t') 
for t — t' — > CP, where 0~ is a negative infinitesimal, is 
related to the population of the level with spin projection 
cr, 



= (4W4(t)) = -«; + (t-OUo- 



(5) 



Therefore our goal is to calculate that particular Green's 
function. As already mentioned, the Hamiltonian in the 
HF approximation is quadratic, so that D~ + (t — t') can 
be calculated either using the re-summation technique 
of Ref.^4[> or the non-equilibrium coherent state func- 
tional integrals [TEl Il6j . We choose the latter method 
because of its relative simplicity. The idea is to define a 
non-equilibrium counterpart of the partition function as 
a functional integral over all fields entering the Hamil- 
tonian and integrate out all degrees of freedom apart of 
those corresponding to d a operators. The resulting ex- 
pression will be quadratic in d CT 's which enables all corre- 
lation functions to be simply read off. Let us define the 



following coherent state functional integral correspond- 
ing to a Keldysh partition function of the subsystem with 
spin orientation <j, 



Z„ = 



J Vd* a Vd a Yl V^* a V^ la exp { - i Si 



i I dt 
lc 



d* a (t)(id t + A + Un_ a )d a (t) 
+ ^7i[<C(*)^(*)+C(*K(*)]]}. (6) 

i 

where by abuse of terminology d* a , d a and ip* a , ipio- stand 
for the coherent state variables of the corresponding oper- 
ators. S{ are the actions of the electrons in each lead. No- 
tice that all time integrations are to be performed along 
the Keldysh contour C. Since the full action is Gaus- 
sian in the ipi a fields we can integrate them out thereby 
generating a non-local action for the impurity electrons. 
Separating time integrations along the C± branches we 
obtain the following expression, 



T>d* a T>d tJ exp ■ 



dtdt' 



x^tJD-^t-O&tO}. (7) 
with d/j denoting a vector of d a (t) for the time variable 
t on C~p D CT (t — t') is the quantity we are looking for 
- it is a matrix components of which give the Keldysh 
Green's functions of the impurity level operators. After 
the Fourier transformation its inverse is given by 



d- 1 H = 



w + A + tfn_ ff + £i1?Gk ( w ) 



ZilfGr+iu) 



(8) 



Therefore the Green's function that we are interested in 
is given by 

D-» = -|D ff Hr 1 ^ 7l 2 G iff +H. (9) 

i 

\D a (u))\ denotes the determinant of the matrix (0. In- 
serting Eqs.(@} into © yields 

(uj + A + U n- a y + (Fi + T 2 ) 

where T; = npo^/f- Finally, according to Eq. JJJ, n a is 
just the energy integral of I|1U|) with an appropriate pre- 
factor, 

8 , /A- V/2 + Un- 
n a = - cot 

7T \ 1 



1-8 



■ cot 



A + V/2 + Un. 



(11) 



where T = Ti + T2 and 8 = T±/T is the asymmetry pa- 
rameter. Eas. i|ll|) are the self-consistency equations for 
the impurity level population probabilities. In the equi- 
librium case these equations reduce to Eqs.(27) of Ref.0]. 
Introducing dimensionless parameters x = —A/U, y = 
U/T and z = V/2T the above equations simplify to 



= — cot 
1-<S 



l [y(r 



x) 



7T 



cot 1 [y{n- a — x) + z] 



(12) 



It is reasonable to expect that both the magnetic {n a ^ 
n_ CT ) and the non-magnetic (n a = n_ CT ) regimes still exist 
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for small z. The critical ix-y) curve, however, acquires 
an additional dimension (z) and thus becomes a critical 
surface. Our goal is to determine its shape. Proceeding 
along the lines of Ref.0, we derive the equations for 



the surface. The first equation is essentially given by 
Eci. (|12|) for n a = n_ CT = n c and the second one is just its 
derivative with respect to n c , 



nric = S cot 1 [y c {n c - x c ) — z] + (1 
k/Vc = 5{1 + [y c (n c - x c ) - z} 2 }- 1 



5) cot 1 [y c {n c - x c ) + z)] , 
{\-5){l + [y c {n c -x c ) + z] 2 }-\ 



(13) 



where the subscript 'c' stands for the parameter of the 
critical surface. Obviously, in the equilibrium situation 
(z = 0) the asymmetry does not matter and one obtains 
the result of Ref . , 



x c = n c — sin(27m c )/27r , 
Tt/y c = sin 2 (7m c ), 



(14) 



where n c plays the role of a parametrisation and runs 
from to 1, see Fig^ Out of equilibrium, the equations 
(|13|l cannot be reduced in the same way and need to be 
solved numerically in the general case. Yet it is possible 
to investigate some interesting special cases analytically. 

Let us first concentrate on the symmetric case 5 — 1/2. 
Then Eas. l|13[) are obviously symmetric with respect to 
transformation z — ► — z, so that we can restrict ourselves 
to positive z. Next we observe that, just as in equilib- 
rium, there is an additional symmetry with respect to the 
axis x = 1/2: x c 1 — x c (which, of course, corresponds 
to n c «-> 1 — n c ). Moreover, at the special point x c = 1/2 
Eas. l|13|l can be solved exactly. The first of them fixes n c 
to 1/2 but does not impose any restrictions on z, so that 
from the second equation we obtain, 



Uc 



7T(1 



(15) 



The main feature of this result is that, at least on the 
line x c = 1/2, the magnetic phase is never completely 
destroyed by the applied voltage. Another remarkable 
property of the phase diagram which can be extracted 
from the result l|15l) is the fact that although the point 
x c = 1/2 on the equilibrium curve corresponds to the 
lowest critical interaction strength U and therefore to 
the most favorable condition for the magnetic phase to 
emerge, it does not possess that property above the crit- 
ical value z* = l/v3- This can be seen from the lowest 
order expansion of Eq. (|15f) in powers of x c — 1/2, 



Uc 



7T(1+Z 2 ) 



l + _(l_ 3z 2 )(a , c _ 1/2) 2 



(16) 



which means that the phase diagram develops a dip to- 
ward the x axis at z > z*, see Fig^ 

Let us now analyse the behaviour of the critical curves 
for small x c . In that limit n c happens to be small as well 
so that we can regard it as a good expansion parameter. 
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FIG. 1: Mean-field phase diagram of the AIM in the sym- 
metric case 5 = 1/2 evaluated numerically using Eqs. ll3H for 
different voltage parameters z. From right to left along the 
Xc — 1/2 line: z — 0, 0.5, 1, 5, 10, 50. The magnetic phase is 
bounded by the curves and the axis of ordinates. 



From the first of Eqs. (|13fl one can easily see that in such a 
situation the cot -1 functions are vanishing so that it can 
be solved after an asymptotic expansion of trigonometric 
functions. The result is y c (n c — x c ) = (7rn c ) _1 + irn c z 2 . 
It can be plugged into the second equation and the whole 
procedure results in the following parametric representa- 
tion, 



x c = n c - [(ir 2 n c ) 1 + n c z 2 ] (7r/y c ) , 

t/z/c = E {! + li^nc)- 1 + irn c z 2 + pz} 2 }- 1 
P =± 



(17) 



Being plotted numerically this curve starts from the ori- 
gin and never crosses the x c = axis again unless 
z > 0.707 w l/v2- Therefore in this latter situation 
the finite voltage encourages the system to become mag- 
netic. This fact also reveals itself as a 'squeezing' of the 
magnetic region toward the lines x c = 0, 1, see FigQ] 

To summarise the effects of the voltage in the sym- 
metric case (5 = 1), the equilibrium conditions for the 
stability of the magnetic phase remain qualitatively in- 
tact as long as z < z*. For larger voltages, z > z*, 
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FIG. 2: Mean-field phase diagram of the AIM in the asym- 
metric case S = 0.4 evaluated numerically using Eqs. ll9> . 
Dashed line corresponds to high negative voltages and dot- 
dashed one represents the limit of high positive voltages. For 
comparison: the solid line is the equilibrium 2 = critical 
curve. 

it becomes favourable to measure the dot level A from 
either of the two shifted chemical potentials [Lb..l in or- 
der to preserve the local moment. Hence the dip in the 
phase diagram. Yet, when V is sent to infinity, the mag- 
netic phase completely disappears in the symmetric case. 
Surpisingly, the magnetic phase persists up to V = oo in 
the asymmetric case (8 ^ 1), which we investigate next. 

The limit of large applied voltage can also be investi- 
gated analytically. However, from the physical point of 
view it is more convenient to interpret the results for a 
system where the voltage is applied asymmetrically. Such 
situation takes place e.g. if we shift our voltage variable 
to zero in one channel and to 2z in the other one. Send- 
ing the voltage to infinity simplifies the first equation of 
the system (TT3)) to 

i \ f n c — 5\ 
y c {n c - x c ) = cot 7T- - , (18) 



from which we immediately see that S < n c < 1. Inser- 
tion of Eq. (fT%|) into the second one of the system lfl"3|) 
yields the following parametric representation, 

( x c = n c - (1 - 5) sin[27r(?v - 5) /(l - 5)}/2n , 

< (19) 
I n/y c = (1-5) sm 2 [n(n c - S)/(l - 6)} . 



At large negative voltages the picture is changed only by 
the transformation x c — > 1 — x c . Sending the applied 
voltage to infinity effectively corresponds to a removal 
the weaker lead. As one would expect in his case the 
critical curve has then exactly the same shape as the 
curve in equilibrium apart of a renormalisation by the 
asymmetry parameter S and a shift of the dot population 
n c , see Fig|5| 

To conclude, we investigated the mean-field phase di- 
agram of the non-equilibrium Anderson impurity model 
by means of the unrestricted Hartree-Fock approxima- 
tion. We present an analytic expression for the critical 
curve between the magnetic and non-magnetic phases in 
the symmetric case when the level energy —A is equal to 
the half of the on-site interaction constant U. It turns 
out that above some critical value of the applied volt- 
age V > 2r/\/3 the critical curve acquires a dip which 
means that the most favorable conditions for the mag- 
netic phase to form change. We show that the magnetic 
phase exists at arbitrary voltages in the asymmetric case. 
We consider our HF solution as an appropriate starting 
point for further investigations into the problem includ- 
ing fluctuations in the functional integral approach or 
through a suitable generalisation of the LMA. 
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